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The relaxation of electrons in quantum dots via phonon emission is hindered by the discrete 
nature of the dot levels (phonon bottleneck) . In order to clarify the issue theoretically we consider a 
system of N discrete fermionic states (dot levels) coupled to an unlimited number of bosonic modes 
with the same energy (dispersionless phonons). In analogy to the Gram-Schmidt orthogonalization 
procedure, we perform a unitary transformation into new bosonic modes. Since only N(N + l)/2 of 
them couple to the fermions, a numerically exact treatment is possible. The formalism is applied to 
a GaAs quantum dot with only two electronic levels. If close to resonance with the phonon energy, 
the electronic transition shows a splitting due to quantum mechanical level repulsion. This is driven 
mainly by one bosonic mode, whereas the other two provide further polaronic renormalizations. 
The numerically exact results for the electron spectral function compare favourably with an analytic 
solution based on degenerate perturbation theory in the basis of shifted oscillator states. In contrast, 
the widely used selfconsistent first-order Born approximation proves insufficient in describing the 
rich spectral features. 



I. INTRODUCTION 

Since the development of quantum well lasers there 
have been continuous attempts to manufacture laser 
structures with even more reduced dimensions. The idea 
behind was to increase the efficiency by enhancing the 
density of states. However, zero-dimensional quantum 
structures (quantum dots) are characterized by a dis- 
crete spectrum, and the recombination probability does 
not depend on the radiative rate alone. Rather, the re- 
laxation pathway into the groundstate becomes decisive. 
Looking at the nearly monoenergetic longitudinal-optical 
(LO) phonons, an efficient relaxation between two dot 
levels seems to be possible only if level distance and LO 
energy match (resonance condition). This type of argu- 
ment has been called phonon bottleneck. P, Q 

^From the experimental side there is an ongoing in- 
tense debate on whether or not the phonon bottleneck is 
seen in the data. 3] However, the recently found ground- 
state lasing in quantum dots under cw conditions seems 
to prove that the phonon bottleneck is not an obstacle 
when trying to increase the laser efficiency by dimen- 
sional reduction. Q 

Nevertheless the theoretical concepts are still contro- 
versial. Obviously, the bottleneck argument relies on the 
assumption of strict energy conservation in the electron- 
phonon scattering, as dictated by Fermi's golden rule. 
A next step towards a realistic description seems to in- 
corporate the intrinsic lifetime broadening of dot levels. 
Krai et al. [f| went along this way by calculating the com- 
plex electron selfenergy due to LO-phonon interaction. 
They claimed that the convolution of initial and final 
state spectral functions gives rise to a broadening which 
is able to circumvent the phonon bottleneck. Arakawa 
and coworkers have treated the electronic transition 
and the LO modes in closed form by wave function evo- 



lution. They pointed out that the final decay of the LO 
phonon into acoustic phonons is decisive for the relax- 
ation process. A different argument includes Auger-like 
electronic excitation for overcoming the sharp ener gy s e- 
lection inherent to the LO-phonon relaxation, pi, |{Olo| 

A more general question concerns the appropri- 
ate theoretical tools for describing relaxation in zero- 
dimensional systems. Non-equilibrium Green's functions 
are often too tedious to be used in realistic models. 
Therefore one is tempted to look for the one-particle 
Green's functions and their broadening as a signature 
for relaxation. 0, ITlj However, standard selfenergy ap- 
proaches as the selfconsistent first-order Born approxi- 
mation p| have to be questioned since they rely implic- 
itly on the existence of an electronic continuum which is 
missing in zerodimensional systems. 

In this paper we want to look closely at this question 
and to qualify the standard approximation schemes (as 
the selfconsistent Born approximation) in application to 
quantum dots. To compare with we present results from 
an exact diagonalization of the electron-phonon Hamil- 
tonian. This can be achieved even for an unlimited num- 
ber of phonon modes provided they have no dispersion, 
and uses a unitary transformation among the phonons. 
If N electronic dot states (fermions) are considered, only 
iV(iV+l)/2 of the new modes (bosons) couple to the elec- 
trons, and for moderate numbers of N the transformed 
Hamiltonian can be easily diagonalized numerically. An 
upper limit of the boson occupation numbers can be fixed 
in accordance with temperature. For N = 1, this exact 
solution is known for a long time as independent Boson 
model, [l^l To the best of our knowledge, the extension 
to a finite number of levels with the important inter- 
level coupling (phonon transition) is presented here for 
the first time. A single boson mode in resonance with 
an equidistant series of electronic levels has been treated 
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exactly in Ref. predicting the phonon staircase effect. 

A related problem is the electron-phonon coupling in 
semiconductor point defects. Both the internal defect 
transitions in the near infrared and the phonon satellites 
of interband transitions show a rich spectrum. Even away 
from strict resonance with an electronic transition, the 
measured transition energies differ from the bare phonon 
value. This has been called LO phonon-donor bound 
state, similar to the exciton-phonon complex introduced 
earlier. [Tij In first attempts for a quantitative under- 
standing, perturbation theory has been used.^^] The 
crystallographic symmetry of the defect dictates what 
kind of lattice distortion (local phonon mode) couples to 
the electronic transitions. Taking into account only a 
few of these symmetry-adapted lattice modes, a full nu- 
merical diagonalization is possible nowadays. |16j This is 
in particular important when dealing with strongly po- 
lar material. We show that semiconductor quantum dots 
behave similarly with respect to the lattice coupling, un- 
derlining once more that the quantum dot behaves as a 
kind of mesoscopic atom. Note, however, the quite dif- 
ferent length scales involved. Whereas the local modes in 
the defect problem are constructed using large parts of 
the Brillouin zone, only a minor part around the T point 
is involved in the quantum dot case. Consequently, the 
LO phonon dispersion can be safely neglected here. 

In Section II the method is outlined and applied to the 
most simple but nontrivial case of two dot levels. Re- 
sults for the spectral function are given in Section III. 
It shows a kind of avoided level crossing if level distance 
and phonon energy nearly coincide, which resembles the 
phonon polariton feature. Still, the exact spectral func- 
tion consists of a series of sharp lines. The selfconsistent 
first order Born approximation fails in this respect by ex- 
hibiting broad spectral features. However, we are able to 
derive a simple analytical approximation which almost 
coincides with the exact results. This employs degener- 
ate perturbation theory for those electron-phonon states 
which are strongly coupled near resonance. It is called 
rotating wave approximation since it sesemble a simi- 
lar treatment of the Jaynes-Cummings model in quan- 
tum optics. Some consequences of the present work on 
the general description of relaxation in quantum dots are 
given as well, and conclusions are drawn in Section IV. 
In the Appendix, the coupling constants are calculated 
adopting parameter values for an idealized GaAs quan- 
tum dot with parabolic confinement. 



II. THE MODEL AND THE 
TRANSFORMATION 

Let us consider N discrete electronic levels i = 
... N— 1 coupled to M phonon modes qi . . . qM of fixed 



energy huiQ. The Hamiltonian reads 

H = Cj cj Cj + hui 6 q bq 

i q 
+ ^^(fiq + feljcjc,-. (1) 

Here, the Ci(c\) denote the fermionic creation (annihi- 
lation) operators, respectively, and & q , 6 q are the cor- 
responding bosonic operators. The coupling constants 
between phonons and electrons, M q J , depend explicitly 
on the fermionic states involved (transitions between dot 
levels). For H to be hermitean, {M^)* = Ml' q must 
hold. For simplicity, the spin is neglected because it is 
conserved by the electron-phonon interaction. 

Since all M bosonic modes couple to the electronic 
levels, a straightforward diagonalization is not feasible. 
We proceed by mapping the phonon operators {6 q } onto 
a new set of bosonic operators {B\}, with the goal that 
only a limited number couples to the electrons. 

We start with an arbitrary linear combination of the 
operators & q written as A\ with A = 1 . . . M, which also 
span the bosonic Hilbert space. Following the well-known 
Gram-Schmidt orthogonalization procedure, 

_ / - - 

B\=A x -^[A x ,Bi a ]B a , B X = B X /^J[B X ,B X ] (2) 

a=l 

we arrive at a new set B\ whose members obey the canon- 
ical Bose commutation relations, [B x ,By] = S\,\>. Then, 
the transformation matrix U with B\ = ^ q U\.qbq is 
unitary, and we have 

E B i^ = E 6 q v ( 3 ) 

A q 

Since the bosonic modes have the same energy, the last 
relation shows that the free boson term in the Hamilto- 
nian Eq.Q remains diagonal. 

The first linear combinations are chosen as 

•1a,.. X^.'.'V ( 4 ) 
q 

where A(i,j) runs over the N(N + l)/2 pairs with 
i > j. The remaining A\ can be taken arbitraryly but 
linearly independent. It follows from the prescribed one- 
to-one mapping that the electron-phonon interaction con- 
tains only the restricted set of the new operators. 
This reduces the numerical labour enormously since the 
relevant Hilbert space contains now N fermionic and only 
N(N + l)/2 bosonic degrees of freedom. 

A. Reduction to a Two-Level System 

The model can be used to describe electrons in a quan- 
tum dot which are coupled to LO-phonons. Under the 
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assumption that the third electronic level is energeti- 
cally well above the lowest two, we will limit the num- 
ber of states to N = 2. Choosing A x = ^ q M^°ti qj 

A 2 = J2q M q°K A 3 = Eq M q 1& q W(3 accomplish that 

the transition matrix element couples only to three non- 
trivial bosonic modes. We obtain 

H = e clc + <l\c\ci + ^TiujoBIBx (5) 

A 

+ (C^+CIB^ (cjco + cjei) 

+ (c 2 Bt + C* 2 B\ + C 3 B 2 + C^l) cJ Co 

+ {CiB x + CIB\ + C 5 B 2 + CIB\ + C 6 B 3 + QB^ c\c x 

with the six coupling constants Ci which follow from 
the transformation J3J). If we further assume that the 
electronic wavefunctions in the quantum dot exhibit a 
well-defined parity the constants C 2 and C4 vanish, and 
B\ couples only to the transition — > 1. The model 
will therefore show prominent features of the Jaynes- 
Cummings model which has been introduced to de- 
scribe a single photon mode coupled to an atomic transi- 
tion. In the present case, however, the detuning has to be 
defined as A = Tiuj — (e\ — eo). Further details regarding 
the explicit calculation of the coupling constants Ct are 
given in Appendix We adopt a parabolic confinement 
potential with extensions yo — zq and xq > yo, having 
in mind an anisotropic harmonic quantum dot. In this 
case, all the lowest dot levels have equal energy separa- 
tion, and the mentioned truncation to just two levels is 
not realistic. Nevertheless it will be applied here to keep 
the numerics at a reasonable level. 

By means of an appropriate rotation of the operators 
B 2 and B 3 , i. e. 



B 2 



((C 5 
(- 



■c;b 2 



C 3 )B 2 h 

(c; 



C e B 3 ) h 

-c;)B 3 )h 



(6) 



with 7 2 = |C 5 - C 3 



I C 6 1 2 , we manage that the new 



mode B 3 only couples to the fermionic particle number 
operator Af = CqCq + c\c\. If we leave out the bosonic 
modes which do not couple to the fermionic levels at all 
we obtain 



H = e cjc + ejcjci (7) 

+ Tiuj B[B x + {CxBx + ClB\){c{c + cjci) 

+ HuoBlBa + (7S2 + iB\)c\ Cl - (r/*5 2 + ijB^Af 

+ hu) BlB 3 - {k*B 3 + kB\)N 



with the new parameters 

r ? = C 3 *(C 3 -C 5 )/ 7 , * = C%Cth 



(8) 



The Hamiltonian conserves the electron number, 
[H, Af] = 0, and consequently the Hilbert space can be 
decomposed according to the electron number (zero, one, 
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FIG. 1: Left: The parameters Ci, 7, 77 and k in units of 
the LO energy as a function of the dot size xo. Right: The 
relative polaronic shifts AE N and the relative shifts of the 
level spacing Ae^. 



or two). In the subspace of zero or two fermions, the 
Hamiltonian can be diagonalized easily because the non- 
diagonal transitions <-» 1 are impossible here. More 
demanding is the subspace of one fermion, which can, 
however, be rationalized a lot by introducing shifted op- 
erators 



B 2 = B 2 -r,N 



ic\ci 



B 3 = B 3 - kM 



(9) 



Note that from Eq.@ onward, the LO energy %uio is 
taken as unit of energy in the remainder of this section. 
We want to stress that the shifted bosonic operators B 2 
and B 3 still obey the canonical commutation rules, but 
do not commute with the fermion operators. The Hamil- 
tonian Eq.Q is now represented as 



H = 



b\b x h 



c[ci — E: 



p 



b\b 2 



b\b 3 



(10) 



ClB\){c\c + cl Cl ) , 



with the modified energies depending of the number of 
particles 



ef = ei + 2Re(777*)AA-| 7 | : 



<N 



l )AT 2 



(11) 
(12) 



The left hand side of Fig. ^ shows the (dimensionlcss) 
parameters calculated for a parabolic confinement poten- 
tial with extensions of yo = zq = 3 nm as a function of Xq . 
All input parameters refer to GaAs, e.g. Huq = 36.7 meV 
(see AppendixEJ. The coupling constants are almost in- 
dependent of the dot size and of order 0.1haj . This value 
of 3.6 meV compares favourably with the polaron shift of 
electrons in bulk GaAs, 2.3 meV. 

The right hand side of Fig. shows the polaronic 
shift due to the phonons B 2 and B 3 , Ea. l|12|) . and the 
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renormalisation of the level spacing Eq. l|ll|l . Since both 
depend on the number of electrons present, we display 
the relative shifts when adding one electron, AE N = 



E 



N 



E 



N-l 



and Ae N = ef - ef 



~N-1 



B. Solution of the Model 

The eigenvectors \N; n 1; n 2 , n 3 ) for the TV-electron 
Hilbert space with N = 0, 2 are simple and given by 

|0; m,n 2 ,n 3 ) = \V) |m)S |n 2 ) 2 |n 3 )g (13) 
|2;ni,n 2 ,n 3 ) ee cj cj | V") |m) J |n 2 ) 2 -2r )+7 l"3)_ 2 « 

with the corresponding eigenvalues 

<,„ 2 ,n 3 = ( £ o + ~ E 2 p )N/2 + n 1+ n 2 + n 3 . (14) 

Here, |V) denotes the electron vacuum, and 



\n)^(n\)-^(Bl + aT\0y c 



(15) 



are the shifted oscillator eigenstates (or coherent states) 
whose vacuum is defined as (B x + a)|0)* = 0. We will 
omit the upper index A = 1,2,3 of the bosonic states 
from now on since confusions are unlikely. 

Considering the N = 1 Hilbert space the Hamilto- 
nian can be easily diagonalized numerically since only 
two bosonic modes are involved. An approximate ana- 
lytical solution is possible as well, where the transition 
matrix element is treated in degenerate perturbation the- 
ory using the coherent states Ea. (jT5)l as basis. In oder 
to show the basic features of this analytic solution we 
neglect for the moment the third bosonic mode B 3 which 
reduces the number of quantum numbers that have to be 
kept track of. Since this mode commutes with the rest of 
the system, it can easily be incorporated afterwards. In 
the same spirit we neglect the polaronic shift l|12l) and the 
shift of the level spacing (jl 1|> which are restored, however, 
in the numerics. 

In the resonant situation (zero detuning A) and ne- 
glecting terms of the order 7 2 Ci, degenerate pertubation 
theory leads to the level repulsion (or avoided level cross- 
ing) 

E ni,n 3 = £ + n l + "2 ± V^l lo(«2|n 2 ) 7 C\\ , (16) 

for ri\ > 1, n 2 > 0. 

The overlap integrals of two displaced oscillators ap- 
pearing in Ea. (|16J) are also known as Franck-Condon fac- 
tors which were first introduced in the theory of excited 
molecules. More generally they are given by 



o(n|m). 



7 



n+m 



-7 2 /2 



\/n\ ml 

min(n,m) 
fc=0 



(17) 



-2k 



kl (n — k)\ (m — k)\ 



and are related to the associated Laguerre polynomials. 
Since the Franck-Condon factors 0(^21^2)7 are oscillat- 
ing as functions of n 2 we find a complex level structure if 
many bosons are present, i.e. at elevated temperatures. 
Introducing a small external broadening will lead to a 
continuous but still highly structured spectral function. 

The energy splitting of Ea. (|16() can be recovered in a 
truncated Hamiltonian where only the nearly resonant 
transitions — > 1 with phonon absorption and 1 — > 
with phonon emission are kept, 



H 



,:W =e clc + e 1 c\c 1 +B\B 1 

+C 1 B 1 c\c + ClB\cla . 



B\B 2 



(18) 



In analogy to the optical equivalent in the Jaynes- 
Cummings model, we use the term rotating wave approx- 
imation (RW). H RW can be represented by 2x2-matrices 
[l8| in the basis 

K,n 2 )° ee cj \V) K) |na>_„ (19) 
\ni,n 2 ) 1 = c\ \V) \n\ - 1) |n 2 )_ r)+7 . 

Notice that the shift of B 2 depends on the electronic level, 
a feature that gives rise to the Franck-Condon factors. 
The corresponding eigenfunctions are 



|ni,n 2 ) ± = c* ina |JM . n 2 
with the eigenvalues 



\ni,n 2 ) 



(20) 



E 



" '2 .1 (^O 



ei) + m + n 2 - 1/2 ± Rn un j2 , (21) 



with the Rabi splitting n2 = A 2 + ni 10(^2^2)7 C\ | 2 
and the weight factors c± i n2 = [(i?„ 1 ,„ 2 ±A)/2i?„ 1 ^ 2 ] 1 / 2 . 

In this approximation, H RW has been mapped onto 
the Jaynes-Cummings model plus an additional bosonic 
mode which merely renormalizes the coupling constant, 
Ci -> {n 2 \n 2 ) 1 C\. 



III. THE SPECTRAL FUNCTION 

In this section we calculate the spectral function of our 
system. We contrast two different approaches and also 
provide the full solution via numerical diagonalization of 
a matrix spanned by only two bosonic modes. 



A. Analytic Expressions 

The spectral function of the electron-phonon Hamil- 
tonian Eq.Ql is usually calculated by means of finite- 
temperature Green's functions (Chapter 6 in Ref. lT^ . 
The selfconsistent first-order Born approximation for the 
(retarded) selfenergy is often used provided the assump- 
tion of weak coupling holds. With some simplification 
it has been applied to the present quantum dot problem 
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with two levels by Krai et ai.jjj There are two diagrams 
of first order in the phonon propagator which can be clas- 
sified as Hartree and exchange selfenergy, E = Tj H + E x . 
Usually, the Hartree term is neglected in view of a con- 
stant electron charge density which is compensated by 
a positive background. However, in the case of quan- 
tum dots being localized in space, the situation is differ- 
ent: In a phonon- assisted transition between levels, the 
charge structure of the electron (given by the confinement 
functions) changes, and a classical electrostatic contribu- 
tion to the lattice deformation appears. Explicitly, the 
Hartree selfenergy of level i is given by [% = 1) 



0J 



(22) 



The electronic occupations numbers rii have to be deter- 
mined via the spectral function, 



(23) 



Under equilibrium conditions, the electrons are dis- 
tributed over energy according to the Fermi function 
f{u>) = (e^ - ^ + l)^ 1 with inverse temperature f3 = 
1/ksT and chemical potential p. Spectral function and 
selfenergy are related as usually via 



AfH = -ilm 



The exchange selfenergy reads 



1 



(24) 



M = X>«m^/ du'Ajo*') (25) 



N{u ) + f{Lj') , N(u ) + l- f(uj') 
j — ujq — uj' + iO 



iO 



Here, N(oj) = [eP^ — I) -1 is the Bosc function and gives 
the phonon occupation. For the selfconsistent first-order 
Born approximation, Eas. l|22l25l) have to be solved in an 
iterative manner until convergency for the total selfen- 
ergy Si(w) is reached. 

This approximate result wil be compared with the 
spectral function based on the numerically exact eigen- 
functions and eigenvalues. As starting point we use the 
general definition 

Aijiw) = - Re J dte luJt Tr (p [a(t), c](0)]+j (26) 

with p as the density matrix. Expanding into exact eigen- 
states we obtain 



x6(e> 



N JV-1 



■U>) 



(27) 



where H\N;u) — e^\N\v) denotes the eigenvalue prob- 
lem in the iV-electron subspace. In equilibrium, p is 



diagonal, (N; v\p\ N; p) = S^.^p^. Due to the well- 
defined parity of the Hamiltonian Ea. (|10|l . the spec- 
tral function has only elements diagonal in the sub- 
level index. Parity means here that the electron-£?i 
subspace decomposes into two subsets defined by the 
states {c f \V)\2ni)o , c\\V)\2m + 1) } and {4|U)|2m + 
l)o, cj|U)|2ni) }, respectively. 

For the spectral function in the rotating wave approx- 
imation, an explicit result can be given using the re- 
sults of Eas. (|2QI21fl . Again the mode B3 is neglected 
which can be easily incorporated afterwards. To com- 
plete the spectrum of H RW we have to include the non- 
interacting states |0, n 2 }° = c' \V) |0} 1«-2) —7? with eigen- 
values eo+J^o-'O) which resemble the groundstate for fixed 
n 2 . 

The spectral function in RW has two contributions, 
A™ = A™' 10 + A™' 21 , which refer to the transition 
between one/zero particles and two/one particles, repec- 
tively. 



A T- W = E (pU+Po^IoKK),! 2 (28) 

ni ,n r 2 =0 

x S(e + (n 2 — n 2 )uj — uj) 

+ E (Pm,n 2 - + Pm,n' 2 )\ C ni,n 2 \ 2 \o(n' 2 \n2}r,\ 2 

ni — l\n2 : n' 2 — 

x (5(e - A/2 - i? nii „ 2 /2 + (n 2 - n' 2 )uo - lu) 

ni — l;n2,n' 2 —0 

x 6(e - A/2 + i? ni ,„ 2 /2 + (n 2 - n 2 )uj - d) , 

wb -ere p^ n2 = °(0, n 2 \p\0, n 2 )°, pl un2t ± = 
± {n ll n 2 \p\n ll n 2 ) ± and p° liM2 = (0;m,n 2 |p|0;ni,n 2 ) 
[see Ea. (II3l) ]. An analogous expression can be given for 
A ,21 (ui). The spectral function ^4q W (o;) satisfies the 
strict sum rules for zeroth and first moment which give 
normalization and average energy, respectively. 

Since the overlap between displaced oscillators does 
not vanish even for different quantum numbers, the dou- 
ble summation over n 2 and n' 2 will account for the satel- 
lites of the spectral function which appear at multiples 
of the phonon frequency loq. The approximation is good 
as long as the system is close to zero detuning, i.e. as 
long as degenerate pertubation theory works. Treating 
additionally the eigenstates in perturbation theory leads 
to marginal improvements only. 

For the construction of the one-particle spectral func- 
tion, a combination of states referring to different fermion 
subspaces are needed, since the spectral function gives 
the frequency-resolved probability for removal (or ad- 
dition) of one electron. A signature of this general 
behaviour is the occurence of both the Rabi splitting 
and the bare detuning in the delta function argument 
of Eq. (|28|l , accompanied with the overlap between non- 
interacting and interacting states. Since, on the other 
hand, the phonon-assisted transition of an electron oc- 
curs exclusively in the subspace of one fermion, we have 
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severe doubts on using a convolution of one-particle spec- 
tral functions for describing relaxation, as done in Ref.|j|. 
We expect this deficiency to be dramatic in particular 
when only a few electronic levels are involved, as in the 
present case. 



B. Numerical Results 

In the following figures we present our numerical re- 
sults. The exact spectral function of the lower level Ao(lj) 
(solid line), obtained from numerical diagonalization, is 
compared with the simplified solution in rotating wave 
approximation A^ w (dotted line), and with the selfcon- 
sistent first-order Born approximation A® (dashed line). 
Both Aq and Aq are convoluted with a Gaussian of 
variance a = hojo/70. 

However, we want to emphasize that the exact spectral 
function (and Af™ (uj) as well) consists of delta func- 
tions only. This was clear from the beginning since a 
finite pertubation cannot change the character of the 
spectrum of the unperturbed system. Since we started 
from dispersionless bosonic modes, the discrete elec- 
tronic spectrum cannot be altered by the finite electron- 
phonon interaction. (19j For the density matrix in Eq. I|26|) 
we assume the grand canonical equilibrium distribution 
p oc exp(— (3(H — fiN")) and fix the chemical potential 
fi halfway in energy between the bare electronic levels. 
This is close to putting just one electron into the dot. 

Figure |3 shows the spectral functions at ksT — Twjq 
(T = 426 K) for the bare detuning zero, i. e. e\ = Tilu$ 
(eo is taken as zero of energy in what follows). This level 
spacing corresponds to a dot size of x$ = 3.9 nm. 

The prominent feature is the splitting of the spec- 
tral function into some kind of dublett which can be 
traced back to the level repulsion in the (dominant) B\- 
channnel. Its value is related to the strength of C\ (3.3 
meV). Looking more closely it becomes apparent, how- 
ever, that the upper and lower structure have not an 
equal weight as one would expect for zero detuning. This 
clearly indicates the shift in level spacing (renormalizcd 
detuning). Further, an overall shift with respect to the 
bare energy eo = is obvious which stems from the po- 
laronic shift of the other modes. All these features can 
be seen also at ksT = Tiljq/4: (T = 106 K) which is more 
relevant for realistic dot spectroscopy (Fig. |3J). The right 
hand panels of Figs. [2and[3]show the first satellite struc- 
tures. They are reduced in weight but have a similar 
appearance as in the main structure. 

In Figure 21 the spectral functions at fcgT = Tiloq are 
compared for two nonzero detunings A = Tiloq — t\ + eo 
which refer to dot sizes of xq — 3.6 nm (left) and xq = 
4.4 nm (right), respectively. 

The spectral function Aq W (uj) shows good agreement 
with the exact solution in all features. In particular, the 
complex structure of lines having different weight and 
position is well reproduced. This means that not only 
the eigenvalues are correctly approximated by the reduc- 
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FIG. 2: The spectral functions A (full), A$ w (dotted) and 
Aq (dashed) at fcsT = Hljo and zero detuning. The energy 
scale is in units of huo, with eo = 0. The right panels show 
the first satellite structure at one phonon energy above (top) 
and below the main level (bottom). 
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FIG. 3: The spectral functions at UbT = and zero 

detuning. Other data as in Fig. 2 



tion of the system into (many) 2x2-matrices, but also the 
eigenvectors. 

Turning to the selfconsistent first-order Born approx- 
imation, Aq(u>), the agreement is not satisfying at all. 
At the higher temperature (Fig. instead of the com- 
plex structure only two broad bands are seen which have 
nearly smeared out the level repulsion. For the lower 
temperature (Fig. [3J|, the Born approximation works 
somewhat better since here the unperturbed but shifted 
ground state energy carries the dominant weight. 
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will be followed in future work. 



0.0 



-0.2 0.0 



FIG. 4: The spectral functions at ksT = Tnoo for negative 
detuning (A = —0.2hu>o, left panel) and positive detuning 
(A = 0.2 huio, right panel). Other data as in Fig. 2 



IV. CONCLUSIONS 

We have presented a model that is suitable to describe 
the electron-phonon interaction in quantum dots and can 
be treated numerically exact. Thus we were able to com- 
pare the exact spectral function with the selfconsistent 
first-order Born approximation. The first observation 
was that the exact spectral function consists of delta 
functions whereas the Born approximation gives a con- 
tinuous spectral function. But even after an artificially 
broadening of the exact spectral function, striking differ- 
ences remain which are due to the unability of the Born 
approach to exhbit all shifts and splittings in detail. We 
would like to mention that without including the Hartree 
selfenergy into the Born calculation, the agreement would 
be even worse. 

As an alternative we propose an analytic solution 
which employs the coherent states for representing the 
electron-phonon Hamiltonian. The rotating wave ap- 
proximation allows to reduce the problem to an (infinite) 
number of decoupled 2x2-matrices. Close to resonance 
(zero detuning), this approximation works remarkably 
good. The appearance of Franck-Condon factors points 
to the great similarity of the electron-phonon interaction 
in molecules and quantum dots. 

Starting from the full electron-LO-phonon Hamilto- 
nian with Frohlich interaction, we were able reduce the 
problem in a quantum dot to a few boson modes only. 
Due to this considerable reduction in numerical labour, 
it was possible to construct the exact eigenstates and en- 
ergies. In the present work, these results have been used 
to calculate the one-particle spectral function. However, 
an extension to two-particle expectation values which are 
relevant for e.g. transition rates is possible. Using the 
density-matrix scheme, relaxation processes can be stud- 
ied within the same reduced Hilbert space, a route which 
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APPENDIX A: DETERMINATION OF 
COUPLING CONSTANTS 

With the transformations of Eq. J2J) we obtain the 
following expressions for the six coupling constants: 



q 
q 
q 

]T M^M^/d + C 5 C 2 /C 3 (Af ) 

q 

q 

J2M^M^-\ Ci f-\C 5 f 
q 

2Re [(2C 2 *C 3 C 4 C 5 * - |C 2 | 2 |C 5 | 2 )/C 2 ] . 



The standard Frohlich coupling for the electron-LO- 
phonon interaction is adopted and applied to the dot 
confinement states, 



\Ci\ 2 

C 2 



\c,\ 2 



c 5 



|C fi l 2 



<J,2J — 

q 



A *ii a2 t e 2 / 1 1 

— — $^ , A = Hu; — 

V 1 ' 2 q q 2e V K oc n 



d 3 r^*(r)e lq - r ^(r) 



(A2) 



For simplicity we consider an anisotropic parabolic po- 
tential as dot confinement, with x as the long axis. The 
two energetically lowest wave functions read 



V'o(r) 
^i(r) 



2Ttx y z j 
-V'o(r) 



s -1/2 



1 / X 



exp - 



yl + _ 

4 \x'o ' vl 4 



(A3) 



where a; > yo,z are the spatial extensions (variances) 
of the groundstate. The relevant level distance is 



£i - £o 



t2 -2 

n x 
2m e 



(A4) 
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The matrix elements in Ea. <|A2|) read 

= cxp -1 (q 2 x xl + q 2 y y 2 + q 2 z zl) 

= i9;E * $°< (A5) 

*i >x - (i - • 

The final integration over q leads to C2 = 0, C4 = on 
account of the well-defined parity of the wave functions. 
The other constants can be reduced to the following (el- 
liptic) integrals, 

42 poo 

In = ^75^ J dtt- n [t(t - a y )(t - a z )}- 1/2 , (A6) 

with OL y = 1 — (yo/xo) 2 < 1 and a z = 1 — (zq/xq) 2 < 1. 
We obtain finally 

Ci = (Zi/2) 1 / 2 , C 3 -(/o) 1/2 

C 5 = (I -h/2)/C 3 (A7) 

C 6 = (J - Ji + 3/2/4 - Cf ) 1/2 . 



For the numerical calculations we choose material con- 
stants of GaAs, i.e. LO phonon energy hujQ = 36.7 meV, 
conduction band mass m e — 0.067m , and dielectric con- 
stants Hoc — 10.7, k — 12 A. According to Ea. lA4t . 
a (long) dot extension of xo — 3.9 nm gives resonance 
between level spacing and LO energy. The effective 
dot extensions in y and z direction are taken equal, 
yo = zq = 3 nm, which renders the elliptic integrals (|A6|) 
to be simple logarithmic functions. 
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